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ABSTRACT 

Recentlv lKnutson et aLl (|2010l ) have demonstrated a correlation between the presence of temperature 
inversions in the atmospheres of hot Jupiters, and the chromospheric activity levels of the host stars. 
Here we show that there is also a correlation, with greater than 99.5% confidence, between the surface 
gravity of hot Jupiters and the activity levels of the host stars, such that high surface gravity planets 
tend be found around high activity stars. We also find a less significant positive correlation between 
planet mass and chromospheric activity, but no significant correlation is seen between planet radius 
and chromospheric activity. We consider the possibility that this may be due to an observational bias 
against detecting lower mass planets around higher activity stars, but conclude that this bias is only 
likely to affect the detection of planets much smaller than those considered here. Finally, we speculate 
on physical origins for the correlation, including the possibility that the effect of stellar insolation 
on planetary radii has been significantly underestimated, that strong UV fiux evaporates planetary 
atmospheres, or that high mass hot Jupiters induce activity in their host stars, but do not find any 
of these hypotheses to be particularly compelling. 

Subject headings: stars: activity — planetary systems — methods: statistical 



1. INTRODUCTION 

With more than 70 transiting exoplanets (TEPs) now 
knowr0, it has been become possible to detect statis- 
tically robust correlations between the parameters of 
TEPs and their host stars, which in turn yields in- 
sights into the processes that are important for deter- 
mining the physical properties of exoplanet systems. 
Several correlations have already been noted, including 
correla tions between the masses and orbital periods of 
TEP s (|Gaudi et al.|[2005tlMazeh et al.|[2005l: [Torres et al.) 
12008), between their surface gravities and orb ital peri- 
ods (jSouthworth et al.ll2007HTorres et al.ll2008[) . between 
the inferred cor e mass of planets and the metallicity o f 
their host stars (|Guillot et al.ll2006H&uFrows et al.|[2007l ). 
between S afronov number and the host star metallicity 
([Torres et a l. 2008), and between the radii of planets and 
their average eq uilibrium temperature and host metallic- 
ity (.Enoch et al.. . 20101 

Verv recently iKnutson et al.l (|2010l hereafter KHIlO) 
have demonstrated a correlation between the emission 
spectra of TEPs and the chromospheric activity levels of 
their host stars, as measured from the strength of the 
emission lines at the Ca II H and K line cores. Plan- 
ets with spectra consistent with noninverted tempera- 
ture models appear to be found around high activity 
stars, while planets with spectra consistent with tem- 
perature inversions are found around low activity stars. 
In demonstrating this correlation KHIlO also published a 
catalogue of log i?^^ values for 39 TEPs. This new, ho- 
mogeneous sample enables statistical studies of the rela- 
tionships between stellar activity and the physical prop- 
erties of TEPs. 

In this paper we use the sample of log R'^^ values from 
KHIlO to investigate correlations between stellar activ- 

^ Harvard-Smithsonian Center for Astrophysies, 60 Garden 
St., Cambridge, MA 02138, USA; jhartmanOefa. harvard.edu; 
^ e.g. http://exopl anets.org| 



ity and other TEP parameters. We find that there is a 
significant correlation between log R'^^ ^'^'^ planet 
surface gravity log (/p . A similar correlation between stel- 
lar activity (as traced by the temporal variation in an 
index related to logi?^^) and the mini mum planetary 
mass Mp sin i was previously noted by iShkolnik et al.l 
(|2005f ) for a sample of 10 RV planets, though the au- 
thors deemed the correlation to be only suggestive. Here 
we demonstrate that the logi?^^-loggp correlation is 
robust with greater than 99.5% confidence. 

The structure of this paper is as follows: in Section [5] 
we describe the data and conduct the statistical analy- 
sis to establish the correlation, in Section [3] we discuss 
a potential observational bias which might lead to this 
correlation, and in Section|3]we speculate on the physical 
origins of this correlation. 

2. DATA AND STATISTICAL ANALYSIS 

Table m gives the log(i?Hj^), and \oggp values adopted 
for planets with both parameters measured, together 
with the sources from the literature for the surface grav- 
ity measurements. In all cases we take log(i?^j^) from 
KHIlO. In general we take the surface gravity of planets 
from studies which calculated it directly from observable 
parameters (the transit duration, depth, and impact pa- 
rameter, together with the RV semiamplitude, eccentric- 
ity, and orbital period; see Southworth ct al. 2007) in a 
Markov-Chain Monte Carlo simulation, or we calculate 
it ourselves from the given observable parameters. 

Figure [T] shows the relation between logi?^^ and 
loggp; the existence of a correlation is readily appar- 
ent. To determine the statistical significance of this cor- 
rela tion we use the Sp earman rank-order correlation test 
fsee lPress et al.iri992f ). finding rs = 0.45. The probabil- 
ity that a random sample of size iV = 39 drawn from an 
uncorrelated population would have either rs >= 0.45 
or rs <= —0.45 is only 0.4%, so the significance of the 
correlation is 99.6%. If we exclude the two hot Neptunes 
GJ 436b and HAT-P-llb, which one might expect to 
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have atmospheric properties that are different from more 
massive planets, the two long period planets HD 80606b 
and HD 17156b, and also exclude planets orbiting stars 
which have temperatures outside the range over which 
IorR'v,^ has been calibrated (4200 K < Teg < 6200 K; 
INoves et al.lfTosl . we find a stronger correlation of rs — 
0.68, with a false alarm probability of 0.032%. The sam- 
ple size in this case is TV = 23. 

We have also searched for correlations between log i?^^; 
and other parameters such as planetary mass, plan- 
etary radius, the radial velocity semi-amplitude, stcl- 
lar mass, planetary deii sity, the Safronov number 9 
(|Hansen fc Barman! 1200 7t) . planetary equilibrium tem- 
perature (assuming zero albedo), average stellar flux in- 
cident on the planet, orbital period, stellar metallicity, 
and stellar effective temperature. Table [1] summarizes 
the strength of each correlation for three separate sam- 
ples: 

1. M > 0.1 Mj, a < 0.1 AU, and 4200 K < Teff < 
6200 K. 

2. 4200 K < Teff < 6200 K 

3. No restrictions 

In addition to the log R'^^-\og gp correlation, positive 
correlations with > 99% confidence are also seen between 
logi?Jjj^ and 9, and between log-R^j^ and pp. All three 
parameters (logf/p, 9 and pp) scale as Mp/Rp {n = 1 
for 9, n = 2 for loggp and n — 3 for pp). The correla- 
tions seen between these parameters and log i?HK most 
likely have the same origin. We focus on loggp here 
both because the correlation is slightly more significant 
for this parameter than it is for 9 or pp, and because for 
TEPs log qp can be determined dire ctly from measurable 
parameters (jSouthworth et al1l2007D , while the other two 
parameters are dependent on stellar models, which could 
conceivably introduce a bias if there is a systematic error 
in the stellar models which depends on activity. 

By checking several different parameters for correla- 
tions with logi?^j^, we have conducted several indepen- 
dent trials, and must therefore increase the false alarm 
probabilities to account for this. Several of the param- 
eters are strongly correlated (loggp, 9 and pp, as well 
as K and Mp, and Teq^p and {F)p), so these are not 
completely independent trials; we estimate that the to- 
tal number of independent trials is between n = 6 and 
n = 12. The corrected false alarm probability for the 
log5p-logi?^K correlation is 1 - (1 - 0.00032)" = 0.19% 
or 0.38% for n = 6 and 12 respectively. 

We note from table [T] that while the planet mass Mp 
shows a positive correlation with \ogR[^^ with ~ 97% 
confidence, there is no significant correlation detected 
between log i?HK ^^'^ ^^"^ planet radius Rp . 

3. POTENTIAL SELECTION EFFECTS 

The observed correlation between log i?HK ^^"^ log 5p 
may potentially be due to observational biases rather 
than physical effects. The most obvious relevant obser- 
vational bias is the relative difficulty of obtaining high- 
precision RV observations for high-activity stars. This 
might introduce a selection whereby lower mass planets 
are not detectable if they are orbiting high activity stars. 
Indeed, as seen in Table [U much of the log R'^^-log gp 
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Fig. 1. — Top: Surface gravity of transiting exoplanets vs. 
the chromospheric activity of the stellar hosts, as measured with 
logi?^j^. The values and data sources are given in Table [2] Filled 
circles show planets with A/ > 0.1 Afj orbiting host stars with 
4200 K < Teff < 6200 K, crosses show planets with M > 0.1 Mj 
orbiting host stars with T^ff < 4200 K or T^g > 6200 K (i.e. 
outside the range over which log R'^^^ is calibrated) , and open 
circles show two hot Neptunes with M < 0.1 Mj, which might 
be expected to have different atmospheric properties than more 
massive planets. We also label the two planets HD 80606b and 
HD 17156b, with semi-major axes a > 0.1 AU, whose properties 
are less likely to be influenced by the stellar flux. These planets 
have relatively high surface gravities most likely as a result of re- 
duced stellar insolation. Bottom: Same as above, here we only 
show planets with A/ > 0.1 Af j, a < 0.1 AU orbiting stars with 
4200 K < Teff < 6200 K. 

correlation can be attributed to a correlation between 
the RV semiamplitude K and logi?^j^, with less contri- 
bution from the planet radius. This is what would be 
expected for a selection effect of this form. Figure [H^top) 
shows the relation between these two parameters. If this 
correlation were due only to an observational bias, how- 
ever, we would expect the correlation to continue to hold 
when low mass and long period planets are also included 
(sample 2). The fact that the correlation is reduced to 
rs = 0.28 with a false alarm probability of 17% in this 
case is evidence that this relation may not be due to a se- 
lection effect. The lowest K planet in this sample, HAT- 
P-1 1, was discovered around one of the most active stars 
in the sample. Moreover, while a selection effect might 
explain the absence of sample 1 planets in the lower right 
corner of Figure [IKtop) or Figure (Hbottom), it does not 
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TABLE 1 

Correlation Between log -R^j^ and Other Parameters 





Sample 1 


Sample 2 


Sample 3 


Parameter 




FAP'' 




FAP'' 




FAP*^ 




0.68 


0.032% 


0.47 


1.5% 


0.45 


0.39% 


e 


0.62 


0.17% 


0.36 


7.0% 


0.41 


0.91% 


pp 


0.66 


0.056% 


0.53 


0.57% 


0.43 


0.66% 


Mp 


0.45 


3.1% 


0.22 


28% 


0.23 


15% 


K 


0.50 


1.5% 


0.28 


17% 


0.28 


8.8% 


Rp 


-0.031 


89% 


-0.053 


80% 


-0.21 


21% 


Ms 


0.19 


39% 


0.15 


48% 


0.30 


6.8% 


Tcn,P 


-0.19 


38% 


-0.16 


42% 


-0.19 


24% 


{Fh 


-0.21 


33% 


-0.18 


38% 


-0.21 


19% 


P 


-0.29 


17% 


-0.23 


26% 


0.0096 


95% 


[Fe/H] 


-0.19 


39% 


-0.18 


38% 


-0.059 


72% 




-0.31 


15% 


-0.32 


11% 


-0.18 


28% 



TABLE 2 

Adopted Values for, the Stellar, Activity Index and 
Planetary Surface Gravity 





Incr B' a 
log «HK 


[cgs] 


xxei. lug (/p 


CoRoT-lb 


-5.312 


3.0266'' 


10 


CoRoT-2b 


-4.331 


3.6157*' 


1 


GJ 436b 


-5.298 


3.1070 


8 


HAT-P-lO/WASP-llb 


-4.823 


3.0800 


16 


HAT-P-llb 


-4.567 


3.0500 


23 


HAT-P-12b 


-5.104 


2.7500 


19 


HAT- P- 13b 


-5.138 


3.1088*' 


21 22 


HAT- P- 14b 


-4.855 


3.6200 


27 


HAT-P-lb 


-4.984 


2.9570 


8 


HAT-P-2b 


-4.780 


4.2260 


29 


HAT-P-3b 


-4.904 


3.3100 


3 


HAT-P-4b 


-5.082 


3.0200 


3 


HAT-P-5b 


-5.061 


3.2190 


3 


HAT-P-6b 


-4.799 


3.1710 


3 


HAT-P-7b 


-5.018 


3.3406*" 


4,5,30 


HAT-P-8b 


-4.985 


3.2300 


18 


HAT-P-9b 


-5.092 


2.9910 


12 


HD 149026b 


-5.030 


3.1320 


17 


HD 17156b 


-5.022 


3.8810 


15 






Q QnOQ** 
o.ouyy 






-4.970 


2.9630 


3 


riU 8UDU6b 


-5.061 


4.U2U2 


28 


IrrLo-lD 


A '7'3Q 
-4. / OO 


o.zzUU 


Q 
O 


irrjb-zD 


-4.949 


O.2980 


3 


irrjo-ob 


-4.549 


3.4250 


1 Q 


lrii/o-4D 


-0.1U4 


Z.OOOU 


lo 


WASP 1 9h 


-5.500 


2.9900 


14 


WASP-13b 


-5.263 


2.8500 


9 


WASP- 14b 


-4.923 


4.0100 


24 


WASP-17b 


-5.331 


2.5600 


26 


WASP-18b 


-5.430 


4.2810 


20 


WASP-19b 


-4.660 


3.1900 


25 


WASP-lb 


-5.114 


3.0100 


3 


WASP-2b 


-5.054 


3.2870 


3 


WASP-3b 


-4.872 


3.4200 


2 


XO-lb 


-4.958 


3.2110 


3 


XO-2b 


-4.988 


3.1680 


3 


XO-3b 


-4.595 


4.2950 


6 


XO-4b 


-5.292 


3.3316*' 
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References. — 1. lAlon so e t"an l|200g); 2._ [Gibson et al.l (200J); 
3. Torres ct al. (2008): 4. Pal ct al. (2008); 5. Welsh ct al. (2010): 6. 
Winn et al. (2008); 7. McCullough ct al. (2008); 8. Southworth 1200^; 
9. Skillen et al. (2009): 10. Gillon et al. (2009); 11. Triau d et ahl 
(2009); 12. Shporer ct al. (2009); 13. Sozzotti et al. |200 l); 14 ^ 
Hebb et al. (2009); 15. Winn et al. (2009): 16. Bakos et al. (200^ 
17. Carter et al. (2009): 18. Latham et al. (2009): 19. Hartman et al] 
([2009,); 20 ^Southworth et al. (2009); 21. Bakos et al. (2009); 22. 
Winn et all (120101): 23 . iBakos et al. (2009); 24. ,Ioshi et al. (200^ 
25. [ Hebb et all iwWt): 2 6. ,A nderson et al., (,2piq'l: 27. .Torres ct iS] 
I 2010[ ): 28. IHebrard et al.M201(i ): 29. IPli'et al.ll201(tl : 30 IWinn et al.l 
1 2009) 

^ log i?HK values are taken from KHIIO. 

*' Calculated from and a combination of a/Ri,^ T14, p* , 6, z, and 
Rp/B.* from the given sourcc(s). 

than anticipated, so that the inflation due to the increase 
in stellar luminosity with time is greater than the grav- 
itational contraction of the planet over time, the result 
would be a positive log gp-log R[i^ correlation. 

Another possibility is that strong stellar UV flux in- 
creases the evaporation of hydrogen from the atmo- 
spheres of hot Jupiters, le ading to higher metallicity, 
more compact planets (e.g. iLecavelier des Etangs et ahl 
2OIOI and references therein). If this were the case, we 
might expect to see a correlation between planet radius 
and stellar activity, and no correlation between planet 
mass and activity. The fact that the opposite effect is 



^ The Spearman nonparametric rank-order correlation coefRcient 
*' False alarm probability. These have not been corrected for the 
total number of independent trials conducted by searching for cor- 
relations between different parameter combinations. 

explain the absence of sample 1 planets in the upper left 
corners of these figures, unless high log gp planets are in- 
trinsically less common than low \oggp planets and low 
log i?HK ^'-'^ Jupiter host stars are intrinsically less com- 
mon than high log host stars. 

As an additional check on whether observational bi- 
ases against selecting planets with low K around high 
activity stars could be responsible for the log/i^j^- 
log gp correlation, we compare the log R'^^-K relation 
for the TEPs from sample 1 to the sample of RV plan- 
ets from the Cal if ornia- Carnegie Planet Search presented 
by IButler et all ()20Q6D . We consider here 84 planets 
from that sample with S index measurements and with 
0.5 < (B — V) < 1.2. We use the relations from 
iNoves et al.l ()1984f ) to calculate logR'^^ from {B - V) 
and S for each of t hese stars. We also plot the jitter and 
residual RMS from IButler et al.l ()2006D . While the resid- 
ual RMS increases towards higher stellar activity levels, 
it is generally well below K — 30 ms~^, and thus well 
below the lowest K systems in sample 1. Moreover, a sig- 
nificant number of RV planets with K < 30 ms~^ have 
been detected around stars with logi?^^: > All this 
suggests that poorer RV precision due to increased stellar 
activity is unlikely to explain the absence of low surface 
gravity (or low K) hot Jupiters around active stars. 

4. DISCUSSION 

Assuming that the observed correlation between 
logi?Hj^ and loggp is not due to an observational bias, it 
is not obvious what physical processes might give rise to 
it. It is well-known that log R'hk is a decreasing function 
of age for FGK stars fe.g. lSoderblom et aLIHQQl ). but the 
log 5p-age relation implied from the log gp-log i?HK ^^1^" 
tion is opposite of what is expected-that planets should 
contract with age, and no t expand with age. For ex- 
ample, by interpolating the lFortnev et all ()2007[ ) models 
while accounting for the increase in stellar luminosity 
over time, we find that a coreless 1.0 Afj planet orbit- 
ing a 1.0 Mq star on a 2.5 day period should decrease 
in radius from 1.22 i?j to 1.13 i?j between 300 Myr and 
4.5 Gyr. Models such as these, however, are known to 
underpredict the radii of many hot Jupiters. If the ef- 
fect of insolation on planetary radii is substantially larger 
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Fig. 2. — Top: RV semiamplitude of TEPs vs. the chromo- 
spheric activity of the stellar hosts, as measured with logijjj^,. 
Symbols are as in Figure [T] Bottom; Same as above, here we only 
show TEPs with M > 0.1 Mj, a < 0.1 AU orbiting stars with 
4200 K < Teff < 6200 K (filled black circles). We also show RV 

tnets orbiting stars with 0.5 < {B — V) < 1.2 fromlBu tler et aTl 
106) (red open circles) together with the residual RMS of an or- 
bital fit to each RV planet (blue crosses), and the inferred stellar 
jitter (excess residual RMS not accounted for by the formal RV 
errors; green open triangles) for each RV planet. While the jitter 
and RMS for the RV planets does increase with stellar activity, 
the levels are still generally well below the K values of the tran- 
siting hot Jupiters. Moreover, the RV planets themselves do not 
show a strong correlation between K and logH^^, and indeed nu- 
merous RV planets have been discovered with K < 30 ms^^ and 



> -5.0. 



observed casts doubt on this hypothesis. Moreover, since 
the stellar activity should decrease with age, this hypoth- 
esis does not explain how planets could re-inflatc when 
the activity is lowered. 



Alternatively, the presence of hot Jupiters may induce 
activity on the host star, either by tidally spinning-up 
the star's convection zone, or via a magnetic star-p lanet 
interaction (see the review bv IShkolnik ei~all 120091) . In 
this case stars with high log i?HK may not necessarily 
be younger than stars with lower logi?yj^. Evidence 
that the presence of a hot Jupiter is correlated with 
increased stellar X-r ay acti vity has bee n presented by 
iKashvap etahl (|2008[ ). while |Pont| (|2009D found that hot 
Jupiter host stars may exhibit excess rotation. Other 
investigations have found evidence of magnetic activ- 
ity variations corre l ated with planet properties (e.g. 
IShkolnik et"alll2005t iLanza et al1l2009[ ). For both tidal 
and magnetic star-planet interactions, the strength of the 
interaction increases with planet mass. If the log i?HK" 
log gp correlation is a by-product of a more fundamental 
log R[i^-Mp correlation, one might wonder why the for- 
mer is detected with higher significance than the latter. 
A possible explanation is that log gp is determined di- 
rectly from measurable parameters while Mp is directly 
proportional to the stellar mass Ms, which in turn is 
dependent on stellar models. As a result loggp is gen- 
erally determined with better precision, and presumably 
with better accuracy, than Mp for TEPs. However, by 
simulating data sets with the observed Afp-log -Rhk 
relation and Mp — Rp correlations, assuming the scat- 
ter about these relations is intrinsic, and assuming the 
observational errors for Mp and loggp are realistic, we 
find that there is only a ^ 1% probability of the FAP 
of log gp-log being less than 0.1% while the FAP of 
Afp-log is greater than 1%. Even if we assume the 
true observational error on Afp is ~ O.SAfp, the prob- 
ability is only ~ 6%. It is therefore unlikely that the 
log i?^j^-Afp relation is driving the log i?^j^-log gp rela- 
tion. 

In summary, we have identified a significant positive 
correlation between stellar activity and planetary sur- 
face gravity. As far as we are aware this correlation is 
unanticipated, and its cause is unclear. 
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the HATNet project, via NASA grants NNG04GN74G, 
NNX08AF23G and SAO IR&D grants. 
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